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We construct a series of 2+1-dimensional models whose quasiparticles obey non-Abelian statistics. 
The adiabatic transport of quasiparticles is described by using a correspondence between the braid 
matrix of the particles and the scattering matrix of 1+1-dimensional field theories. We discuss in 
depth lattice and continuum models whose braiding is that of 5*0(3) Chern-Simons gauge theory, 
including the simplest type of non-Abelian statistics, involving just one type of quasiparticle. The 
ground-state wave function of an 50(3) model is related to a loop description of the classical two- 
dimensional Potts model. We discuss the transition from a topological phase to a conventionally- 
ordered phase, showing in some cases there is a quantum critical point. 



I. INTRODUCTION 

Understanding phases with topological order has be- 
come an important theme in condensed-matter physics. 
Well-understood examples of topological fluids include 
the fractional quantum Hall states which arise for elec- 
trons in two dimensions moving in large magnetic fields. 
The existence of quasiparticles or quasiholes with frac- 
tional statistics is a central and striking prediction of the 
theory of the fractional quantum Hall effect and follows 
directly from the nature of the electronic correlations in 
this quantum fluidi^. In spite of its profound concep- 
tual importance, only this week has there been a report 
of experimental confirmation of this startling prediction—. 

Non-Abelian fractional statistics are a fascinating 
property of some fractional quantum Hall states^. Here, 
the wave function depends not only on which particles are 
braided, but on the order in which they are braided: the 
states carry a non-Abelian (matrix) representation of the 
braid group. One of the motivations for the current con- 
sideration of non-Abelian phases is the realization that a 
physical system in a non-Abelian topological phase be- 
haves effectively as a universal quantum computer^Si 

The topological quantum fluids arising in the fractional 
quantum Hall effect have an effective hydrodynamical de- 
scription in terms of Chern-Simons gauge theorie s 8 ' 9 ' 10 . 
Pure Chern-Simons is a topological field theory, meaning 
that its correlators are independent of the position of the 
operators and depend only on topological invariants^. It 
has a vanishing Hamiltonian; the only non-trivial proper- 
ties arise from the braiding of its Wilson and its Polyakov 
loops. In this paper we are mainly interested in theories 
whose ground state is topological, but whose gapped ex- 
citations have non-Abelian statistics. 

Fractional quantum Hall fluids do not have time- 
reversal symmetry, but topological order occurs in mod- 
els with unbroken time-reversal invariance as well. These 
"spin-liquid" phases were originally speculated to be re- 
sponsible for the unusual behavior of the "normal state" 
of high temperature superconductorsiSai 3 " although, in 
spite of much effort both theoretical and experimental, 



there is yet no solid evidence for any spin-liquid state. 
Nevertheless, as a consequence of much (theoretical) ef- 
fort, we know that topological order occurs in the ground 
states in certain gapped systems with time-reversal in- 
variance and reasonably-local interactions^, for example 
in quantum dimer models on non-bipartite latticesi^*i^. 
Morally, these topological phases are equivalent (in the 
sense of asymptotic low energy theories) to deconfined 
phases of an effective gauge theory. The excitations of 
these time-reversal invariant phases do exhibit electron 
fractionalizationiSiiLiSiiS, but the statistics is Abelian. 

Non-Abelian topological phases are even harder to 
come by. Such phases with broken time-reversal sym- 
metry do occur in the fermionic Pfaffian (a.k.a. Moore- 
Read) wave function^ for the v = 5/2 fractional quantum 
Hall effect, and in the v = 1 bosonic Pfaffian stato 20 ' 21 ' 22 . 
The former occurs in models of p-wave superconductors 
as well2£. Field theories with such non-abelian statis- 
tics also have been found 2 ^ 2 ^* 2 ^. In a number of these 
cases, the long-distance physics can be described by a 
Chern-Simons gauge theory, which breaks time-reversal 
symmetry. Non-abelian topological phases, however, oc- 
cur in time-reversal-invariant systems as wellSi The 
resulting effective Chern-Simons theory is doubled to 
restore the time-reversal symmetry. To give the the- 
ory a gap while keeping the topological theory as its 
ground state, one can include the electric-field part of 
the Maxwell term J d 2 x E ■ E in the Hamiltonian 27 ' 28 ' 29 . 
Hence in these topological phases, the ground-state wave 
function is a superposition of configurations of Wilson 
loops in two-dimensional space, while the world-lines of 
the excitations correspond to Polyakov loops in the 2+1- 
dimensional field theory 7 ". For exampl e, the lattice mod- 
els discussed in detail in Refs. 30 31 have a continuum 
description in terms of doubled SU(2)k Chern-Simons 
gauge theory. The resulting configuration space is natu- 
rally associated with the Temperley-Lieb algebra 32 . 

A natural description of the configuration space 
of models in a topological phase is in terms of 
loopo 5 ' 6 ' 7 ' 30 ' 31 . This holds in both Abelian and non- 
Abelian cases. Precisely, each basis state in the Hilbert 
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space of the quantum theory is a loop configuration 
in two dimensions. Quantum dimer models and their 
generalizations^*^ can also be viewed as quantum loop 
gases. 

In this paper we reexamine the problem of non-Abelian 
topological phases by starting with the statistics we wish 
to have, and working backward to construct a model 
which exhibits them. We thus first give an algebraic 
way of characterizing the braiding in both SU(2)k and 
SO(3)k Chern-Simons theories. We show that such a 
braid matrix of a 2+1-dimensional theory is a limit of the 
S-matrix of an associated relativistic 1+1 dimensional 
model, and give an intuitive argument as to why this 
is so. We then show how to explicitly construct quan- 
tum two-dimensional models with these braid relations 
by utilizing the structure of the factorizable 5-matrices 
of integrable 1+1-dimensional models. 

Specifically, we embed the 1+1-dimcnsional model in 
two-dimensional Euclidean space, and find a Rokhsar- 
Kivelson-type quantum Hamiltonian 14 acting on this 
two-dimensional space whose ground state has the prop- 
erties expected of a model with non-Abelian statistics. 
In both cases we discuss in detail, the Hilbert space is 
that of a loop gas: in the SU(2)k case, the loops are self- 
avoiding and mutually- avoiding 30 , while in the SO(3)k 
case, the loops intersect. The latter are thus more akin 
to nets than loops 3 ^. Both these loop gases are associated 
with well-known two-dimensional classical statistical me- 
chanical models: in the SU (2)k case, this is the known as 
the 0(n) lattice model with n = 2cos(7r/(fc + 2)), while 
in the SO(3)k case, this is the Q-state Potts model with 
Q = 4cos 2 (7r/(fc + 2)). The loop expansion of the for- 
mer is well knownS, but the one we utilize for the Potts 
model does not seem to have been discussed before. 

Having an explicit lattice construction of the states 
enables us to construct (reasonably) local quantum lat- 
tice models with these ground state wave functions. By 
studying the statistical properties of the absolute value 
squared of these wave functions, we can investigate the 
correlations described by these quantum states, and de- 
termine if they describe quantum critical points or mas- 
sive (topological) phases. Both here and in the SU{2) 
case 3 ^, the result depends on the level k. 

The paper is organized as follows. In Section [n] we de- 
scribe the algebraic approach to non-Abelian statistics in 
both the SU(2)k case and SO(3)k cases. In Section ITTT1 
we discuss quantum loop gases and their relation to the 
S matrix of 1+1-dimensional integrable field theories. In 
Section IIVI we give an explicit construction of these S- 
matrices for SO(3) and SU(2) braiding. In Section [V] 
we discuss the corresponding 1+1-dimensional field the- 
ories. In Section IVII we discuss lattice models whose 
ground states are precisely the loop wave functions, both 
SU(2)k and SO(3)k, with the desired braiding properties. 
We give a set of specific criteria that 2+1-dimensional 
quantum Hamiltonian ought to satisfy and give an ex- 
plicit construction. In Section fVIII we discuss under what 
circumstances these wave functions are topological and 



when do they describe quantum critical systems. In Ap- 
pendix [3] we give a summary of the Landau-Ginzburg 
description of the 1+1-dimcnsional theories whose S ma- 
trix we use. 



II. BRAIDS AND ALGEBRAS 

Particle statistics, of course, are the effect on the wave 
function when particles are adiabatically transported 
around each other a large distance away. This picture fol- 
lows from the familiar concept of adiabatic particle trans- 
port, developed in detail in the context of the Laughlin 
states of the fractional quantum Hall effect where it fol- 
lows from the Berry phase accumulated during an adia- 
batic evolution of the state with two quasiparticles^. 

Adiabatic particle transport can be represented picto- 
rially by drawing the world lines of the particles, which 
are the paths they trace out in space-time. Since the 
particles stay far apart, we need only study paths which 
do not cross. (Non-trivial braid statistics always require 
the assumption that particles have a hard-core short dis- 
tance repulsion). The world lines of the particles there- 
fore braid around each other. Formally, the set of all pos- 
sible braidings is a group, acting on the space of states 
of the system. Different types of statistics correspond to 
different representations of the braid group. In this paper 
we consider two spatial dimensions, where both Abelian 
and non-Abelian statistics are possible. A system whose 
quasiparticles are associated with a non-Abelian repre- 
sentation of the braid group has a degenerate set of states 
of quasiparticles at fixed positions x%, X2, ■ ■ ■ , Xn- Call 
this space of states V(N). The states in V(N) are lo- 
cally indistinguishable but differ topologically. When a 
quasiparticle is taken around another, states in this de- 
generate subspace are rotated into each other. Since this 
is quantum mechanics, it can take a state to a linear com- 
bination of other states: adiabatic particle transport can 
entangle the states. 

Although we will not require its use, it is useful to note 
that non-Abelian statistics can be implemented in a local 
field theory by assigning to each quasiparticle a charge 
and a flux under a non-Abelian gauge group. Within 
this picture, a quasiparticle carries a representation V 
(i.e. a 'charge') of the gauge group G and a flux g e G. 
When quasiparticle 1 is taken around quasiparticle 2, its 
internal state ipi S V% is rotated to 52^1 by the action 
of the flux i?2 S G which is associated with quasiparticle 
2. It is worth keeping in mind that this picture is not 
gauge invariant, and there is, in fact, no local degree of 
freedom associated with each quasiparticle (since a gauge 
transformation can change ip € V into any other ip' £ V 
and can change the flux g into any other element of its 
conjugacy class hgh^ 1 ). 

Studying the statistics of a 2+1-dimcnsional system 
can effectively be reduced to a two-dimensional prob- 
lem. We project the world lines onto the plane (ignor- 
ing boundary conditions), and call them strands. When 



there are N particles, we have N strands. A braiding 
in 2+1 dimensions results in the crossing of two strands 
in the two-dimensional picture. In this projection, there 
are overcrossings and undercrossings. As long as we are 
only interested in the statistics of the particles, the other 
details of the projection are not particularly important: 
we can move the strands around at will as long as we do 
not remove crossings or create new ones. 

It is useful to think of this collection of strands in 
the plane in a 1+1-dimcnsional fashion. The degenerate 
states of the 2+1-dimensional system correspond to a set 
of degenerate multi-particle states in a one-dimensional 
quantum system. In this 1+1-dimensional quantum sys- 
tem, each strand is the world line of a real local degree of 
freedom. Heuristically, this is the "gauge-fixed" version 
of the model. Consider a configuration at t = — oo (i.e. 
before any of the particles have been braided) where all 
the particles are very far from each other (i.e. all at spa- 
tial infinity). We can thus consider these particles to all 
be on a circle. This circle is our one-dimensional space. 
We can construct the full configuration by a sequence of 
adiabatic braidings as the particles move inward. Since 
we are free to move around the strands as long as we do 
not add or remove crossings, we can then take all the 
particles at t = oo to lie on a circle as well. Thus we can 
project all the original 2+1-dimensional world-lines onto 
a two-dimensional annulus. We can then view the an- 
gular direction of the annulus as one-dimensional space, 
and the radial direction as (Euclidean) time. Thus, each 
configuration in the plane can be regarded as an adi- 
abatic evolution in an equivalent 1+1-dimensional Eu- 
clidean field theory in which each strand (or particle) be- 
longs to a given Hilbert space associated with the species. 

To make this more precise, let us consider the TV- 
particle space of states V(N) on which the braid group 
acts. For now we let this space be the tensor prod- 
uct of N copies of the single-particle space of states: 
V(N) = V® N = Vi ® V 2 ® •• • ® V N . (Later, we will 
see that the actual space of states is a subspace of V® N .) 
In the 1+1 dimensional picture, we can think of V® as 
the space of particles on a circle. The elements of the 
braid group corresponding to overcrossings and under- 
crossings are denoted Bi and BJ~ and are displayed in 
Fig. ^ The subscript i means that Bi is describing the 
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FIG. 1: braid, reverse braid, and Temperley-Lieb generator 

crossing between the ith and the (i + l)th particles, and 
so acts non-trivially on the space Vi ® Vi + i, and with the 
identity on the other spaces Vj, with j ^ 1. For ex- 
ample, the braids in Fig. |2]are described algebraically as 
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satisfy the relations 

BiB i+ iBi = Bi + iBiB i+ i, 

B t Bj = BjBi \i~j\>2. (2.1) 

These relate configurations which are topologically iden- 
tical, as can easily be seen from Fig. |3J 




FIG. 3: Consistency relation for braiding 



If the matrices Bi are diagonal, then the statistics are 
Abelian. For bosons the Bi matrices are all the iden- 
tity; for anyons their entries are phases. In this paper 
we are interested in non- Abelian representations of the 
braid group, so that particles obey non-Abelian statis- 
tics: the wave function changes form depending on the 
order in which the particles are braided. One can give 
explicit matrix representations of the braid group. How- 
ever, it is usually much convenient to study the algebra 
of the matrices involved. In the cases of interest here, the 
statistics of the particles can be obtained directly from 
the algebraic relations the matrices obey, without need 
for their explicit representation. 

A. The SU(2) theory 

A famous one-parameter set of non-Abelian repre- 
sentations of the braid group arises from utilizing the 
Temperley-Lieb algebra^. These representations give 
rise to the Jones polynomial in knot theory2§i2L2&, and 
correspond to the braiding of Wilson and Polyakov loops 
in SU (2) Chern-Simons theory. 

The Temperley-Lieb algebra originally arose as a way 
of relating the Potts models to the six- vertex model. The 
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transfer matrices of these two models (and a number of 
other models) can be written in terms of different rep- 
resentations of this algebra. This means that any prop- 
erties of the model which can be computed from purely 
algebraic considerations will be the same for any such 
model. A generator of the Temperley-Lieb algebra ej acts 
non-trivially on the ith and (i + l)th particles; a useful 
pictorial representation is given in Fig. ^ The algebra 
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de u 



(\j~i\>2). 



(2.2) 



where d is a parameter. These algebraic relations are 
drawn in Fig. ^ From the picture, one can see that d 




= d 
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FIG. 4: The Temperley-Lieb algebra 



These braid relations are those of Wilson loops in 
5 l C/(2)fc Chern-Simons theory whenAi 



d = 2 cos 



k + 2 



or equivalently q — e n ^ k+2 \ The integer k is the coeffi- 
cient of the Chern-Simons term in the gauge theory, and 
is known as the level. It must be an integer to ensure 
gauge invariance of the Chern-Simons gauge theorjiiS. 

The configuration space of Wilson loops on the plane 
is a very large space and arbitrary superpositions of these 
states do not obviously describe a topological ground 
state. It turns out that in order to enforce the condi- 
tion that states be purely topological, these ground states 
must satisfy additional properties which can be enforced 
by means of a suitable projection operator. This is known 
as the Jones- Wenzl projector, which acts on k+ 1 strands 
without crossings. For details on its application in this 
context, see Refs. 00 We will give an algebraic descrip- 
tion of this below. By construction, states obtained by 
means of this projector are topological and do not sup- 
port local low energy degrees of freedom. Conversely, 
unprojected states can describe low energy, even mass- 
less, degrees of freedom, and are unphysical states in a 
topological gauge theory such as Chern-Simons. How- 
ever, unprojected states may describe degrees of freedom 
associated with quantum critical points. 



can be thought of as the weight of a closed loop. 

Representations of the braid group can be found from 
representations of the Temperley-Lieb algebra by letting 

B i = I-qe i , (2.3) 
where / is the identity. This is illustrated in Fig. [SJ The 




FIG. 5: The braid in terms of the Temperley-Lieb generator 



Bi defined in this fashion obey the braid-group relation 
(|2.1|) when 

d = q + q~ 1 . 

It is also easy to check that 

Br 1 =1- q-^i. 

Note that in writing the braid group in this fashion, we 
have resolved the crossings B and B^ 1 in terms of strands 
which no longer intersect. 



B. The SO (3) theory 

Here we discuss another one-parameter set of non- 
Abelian representations of the braid group. These de- 
scribe the braiding in 5*0(3) Chern-Simons theory, in- 
stead of SU(2). The algebras are of course the same; 
the key distinction is that Wilson and Polyakov loops in 
S*0(3) Chern-Simons occur in only integer-spin represen- 
tations. The corresponding representation of the braid 
group is given in terms of the £0(3) Birman-Murakami- 
Wenzl (BMW) algebra^, defined below. This algebra 
has two non-trivial generators Xj and Ej acting on ad- 
jacent strands. 

The SO (3) braid relations can be found by "fusing" to- 
gether two strands obeying the Temperley-Lieb algebra: 
the SO (3) BMW generators Xj and Ej can be written in 
terms of the Temperley-Lieb generators ej. Heuristically, 
the idea is to exploit the fact that a spin-1 representa- 
tion of SO(3) can be found from the tensor product of 
two spin-1/2 representations of SU(2). This statement is 
still true in the "quantum-group" algebra Ug^sl^), which 
is a one-parameter deformation of the ordinary Lie alge- 
bra s?2- One can define an action of U q (sl2) on the space 
of states V ® N which commutes with the ; see Ref. 
for an extensive discussion of quantum groups in the con- 
text of non- Abelian statistics (q there is q 2 here) . In par- 
ticular, to relate the two algebras, first note that e^/c? 
is a projector. In U q {sl2) language, this projects onto 
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the trivial spin-0 representation. The projector onto the 
spin-1 representation is therefore 



Pi = I-\ei 
a 



(2.4) 



so that PiCi = 0. The single-particle space of states Wj 
in the so(3) BMW algebra is comprised of two "fused" 
Temperley-Lieb strands, projected onto the spin-1 rep- 
resentation. In an equation, Wj — P 2 j-i[V 2 j-i ® V 2 j}- 
Pictorially, just think of each strand in the so(3) theory 
as the left-hand-side of Fig. El 
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FIG. 6: Projecting onto the spin-1 representation 



With this identification, the so(3) BMW algebra fol- 
lows from the Temperley-Lieb algebra. Since lines never 
cross in the latter, they cannot cross in the former either. 
When two of the fused strands come near each other, 
there are now three possibilities for what happens, which 
we display in Fig. From the pictures, we read off that 
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FIG. 7: The generators of the 50(3) BMW algebra 



Pij-\ Pij+i e 2 j e 2 j-i e 2j+ i e 2j P 2 j-i P 



Xj — d P 2 j-i P 2 j+i e 2 j P 2 j-i P 2 j+i- 



(2.5) 



These generators act on the two-particle states in Wj ® 
Wj+i. It is straightforward to verify using Ij2.2|l that they 
obey the o(3) BMW algebra. We have 



(E 2 ) 2 = {Q-l)E u 
{Xi) 2 = {Q-2)Xi- 
EiXi = XiEi = (Q — l)Ei. 



Ei 



(2.6) 



where the parameter Q = d 2 . 

Relations involving generators on adjacent sites (e.g. 
EiE i+ iEi = Ei) are straightforward to work out us- 
ing the Temperley-Lieb algebra; they can be found for 
example in Ref. ]4£ Most become fairly obvious after 



drawing the appropriate picture. The relations involv- 
ing only the Ei are those of the Temperley-Lieb alge- 
bra (|2.2(l . but with d replaced here by Q — 1, so that 
closed isolated loops of "spin-1" particles get a weight 
Q-l = d 2 -l = l + q 2 + q- 2 . This factor of d 2 ~ 1 
is obvious from the pictures; the d 2 comes from the d 
from each loop, and 1 must be subtracted because of the 
projection onto spin-1 states (projecting out the singlet). 

The reason we have done all this work is to give us an- 
other representation of the braid group. Namely, defining 



B 



SO(3) 



q 2 X - Xj 



l E< 



(2.7) 



it is straightforward to check using the so(3) BMW rela- 
tions that the Bj satisfy the braid-group relations (|2.1|l . 
X is the identity on the projected Hilbert space W® N \ on 
V® N , we have 1 = P 2 j_iP 2 j +1 (see Fig[7|). One can also 
check that 



(B 



SO(3)x_i 



= q- 2 l - Xj 



q 2 E 3 



Particles with braiding given by £? ■ arise from SO(3) 
Chern-Simons theory. This follows from our construc- 
tion: we basically have restricted the particles to be as- 
sociated with integer-spin representations of U q (sl 2 ); this 
is precisely what one does to go from SU(2) to SO(3). 



C. The Jones- Wenzl projector 

The Jones- Wenzl projector is simply expressed in 
terms the projector onto the representation of spin s 



,([fe+l]/2) 



for all 



of U q (sl 2 ): for a given k it is simply 
j. This projector involves k + 1 strands, so this amounts 
to being able to replace the identity acting on k + 1 
strands with a linear combination of other Temperley- 
Lieb or BMW elements. The necessity of this projec- 
tion is also apparent from the representation theory of 
U q (sl 2 ): when q k+2 = —1, the spin-(/c + l)/2 represen- 
tation is reducible but is indecomposable (it cannot be 
written as a direct sum of irreducible representations). 
Performing the projection avoids all sorts of complica- 
tions such as zero-norm states. 

We have already seen one example of such a projector. 
The projector Pj = I—ei/d is the projector onto the spin- 

1 representation of the quantum-group algebra, so — 
Pj. When k = 1 (d = 1), the Jones- Wenzl projector is 
simply Pj : any spin-1 combination of strands is projected 
out. Therefore the SO(3) theory at k = 1 is trivial. 

A case of great interest is SO(3)3, the "Lee- Yang" 
model. This is the simplest model of non-abelian statis- 
tics, because there is only one type of non-trivial braid- 
ing. The naming arises because the braiding relation for 
the particle in this model is the same as the fusion rule 
associated with the Lee- Yang conformal field theory^ 3 .. 
We have for any k in both the SU(2) and SO(3) models 

1 „ 1 



,(2) 



P 2l -lP 



2j+l 



d? 



(d 2 -2)(d 2 



~) E > 

(2.8) 
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When k — 3, the Jones- Wenzl projector sets = 0. 
In the SU(2) 3 theory, this is a relation involving four 
strands, while in the 5*0(3)3 theory, this relates two of 
the fused strands. This means that in £0(3)3, imposing 
the Jones- Wenzl projector allows us to replace any ap- 
pearance of a generator Xj in favor of X and the Ej : we 
have 

X, = (Q - 2)(T+ Ej) for 50(3) 3 . 

Plugging this back into the braid relation (|2.7[1 and re- 
membering that ~ 2 = q 2 + q~ 2 gives 

Bf 0(3)3 = -q 2 l-q- 2 E,. (2.9) 

Note that for this value k = 3 (and this value alone), 
the generators Ej obey the same Temperley-Lieb algebra 
as the ei, because d = 2cos(7r/5) = d 2 — 1. Moreover, 
<7 5 = — 1 here, so the braid generator (|2.9|l is equivalent to 
that in 12.3(1 . We thus are led to an intriguing result: the 
0(3) 3 theory is almost identical to that of SU (2)3. There 
is one important difference: in 50(3)3, we have already 
imposed the Jones- Wenzl projector, while in 5C/(2) 3 this 
still needs to be imposed. Thus (locally) 50(3)3 with the 
Jones- Wenzl projection imposed is equivalent to SU(2) 3 
without it imposed. 



III. QUANTUM LOOP GASES AND THE S 
MATRIX 

In the preceding section we discussed some of the mar- 
velous properties of particles with non-Abelian statistics. 
Now we discuss how to associate the 5 matrix of a rela- 
tivistic 1+1-dimensional field theory to the type of braid- 
ing discussed above. We argue that knowing the 5 matrix 
in this 1+1-dimcnsional theory allows us to find a quan- 
tum loop gas in two space dimensions where the quasi- 
particles should have these braid relations. 

A natural place to look for models with non-Abelian 
statistics is in quantum loop gases. The reason is that 
if one projects the world lines onto the (spatial) plane, 
one obtains loops: strands cannot end. It is also natural 
from the point-of-view of field theory: in pure Chern- 
Simons theory the only gauge-invariant degrees of free- 
dom are loops. In a 2+1-dimcnsional picture, this means 
it is a good idea to look for a system where the low- 
energy degrees of freedom are loops in the plane, a quan- 
tum loop gas. In a number of cases it has been argued 
that quantum loop gases turn into gauge theories with 
Chern-Simons theories in the continuumSi22^4. The 
excitations can be non-Abelian in a topological phase, 
where the ground state contains a superposition of Wil- 
son loops (loops in the spatial plane) . The excited states 
are Polyakov loops, loops which extend in the time di- 
rection. The quasiparticles have non-Abelian statistics 
when the gauge group is non-Abelian with level k > 1, 
Ref. [U 



To understand these loop gases, it is best to first fo- 
cus on the properties of the ground state. The types of 
ground states we are interested in are liquid states, where 
all local order parameters have vanishing expectation val- 
ues. Such a ground state is a superposition of different 
loop configurations. The basis states \s) of the Hilbert 
space can all be described by some configuration of loops 
in two dimensions. In these models, the wave function \t 
of this ground state can be written in the form 

e -S(s) 

(«l*> = -3- (3.1) 

where S(s) turns out to be the action of the classical 
two-dimensional loop model for the configuration corre- 
sponding to s. Z is the usual two-dimensional partition 
function with weight |(s| , I')| 2 , which is the functional in- 
tegral over all configurations s with weight e^ s ^ s ^ s W. 
However, just because the loops can be non-local does not 
mean the action needs to have long-range interactions. 
Recall, for example, that one can describe all the con- 
figurations in the classical two-dimensional Ising model 
in terms of closed loops of arbitrary length, the domain 
walls. Nevertheless, the interactions are still local. 

Since we are identifying the loops with the quasipar- 
ticle world lines, we need to find a quantum loop gas 
whose ground-state wave function satisfies the appropri- 
ate braiding properties. To make this notion precise, let 
us examine the classical loop gas (i.e. the one with action 
iS(s)) corresponding to the ground state. Now view this 
two-dimensional loop g 1+1-dimcnsional quantum 

system. The loops then can be thought of as world lines 
of particles in this 1+1-dimcnsional system. Their wave 
function is then a vector in the space V® N , just like be- 
fore. The braid generators act in the same way as well. 

To summarize the arguments so far: we project the 
world lines of the particles in 2+1 dimensions down onto 
the plane, so that they form loops. A two-dimensional 
quantum system possessing such particles is a loop gas, 
where the degrees of freedom are loops in the plane. The 
ground-state wave function of the quantum system is 
the expressed in terms of the action S(s) of the corre- 
sponding classical two-dimensional loop gas. Finally, we 
then identify these loops as the world lines of particles in 
the corresponding 1+1-dimcnsional problem. The upshot 
is that by restricting ourself to considering the ground 
state, we have reduced a 2+1-dimensional problem to a 
1+1-dimensional one. Theories for which this construc- 
tion holds are inherently holographic in that the degrees 
of freedom can be naturally projected to a boundary. 

In the 1+1-dimensional theory, when two particle 
world lines cross, the 5 matrix plays the role of the 
braid matrix. In other words, it describes what hap- 
pens to the wave function when the paths of two par- 
ticles cross. Consider the wave function describing two 
particles of momentum and position pi,Xi and Pi+i, 
respectively. The 5 matrix is a matching condition on 
the wave functions for X4 <C Xi+i and X4 3> CEj+i. As 
before, the wave function is a vector in V® N . The two- 
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particle 5 matrix for scattering particle i from particle 
i + 1 acts non-trivially in Vi ® V i+ i. 

^V^Vi+x^i > x i+l) = SiiPiiPi+lWviQVi+xfa < ^+l)- 

Our theories are rotationally invariant in two- 
dimensional space, so the corresponding 1+1-dimensional 
theory is Lorentz invariant. This means that the 5 
matrix depends only on the relative rapidity 9: defin- 
ing pi = msinh^i and Pi+i = msinh# i+1 , we have 
9 = 9\ — 02- We note that the 5* matrix here should not 
be confused with what is usually called the modular 5 
matrix, which governs the braidings in Chern-Simons 
theory and in conformal field theory. 

This correspondence between the braid group and the 
5 matrix has long been known, in the context of knot 
theory^. Representations of the braid group (and the 
resulting knot invariants) can be found by taking a spe- 
cial limit of solutions of the Yang-Baxter equation^. In 
physics, the Yang-Baxter equation arises in integrable 
lattice models and field theories. In integrable lat- 
tice models, the Boltzmann weights must satisfy Yang- 
Baxter; in the former, the 5 matrices of the particles do! 
Note, moreover, that the braid matrices and the 5 matri- 
ces are acting in the same space V® N . So our arguments 
indicate that the braid matrices of the 2+1- dimensional 
theory are a limit of the S matrices of the corresponding 
1+1-dimensional theory. 

It is not difficult to find in which limit this holds. The 
Yang-Baxter equation for the 5 matrix arises from re- 
quiring that the three-body 5 matrix factorizes into a 
product of two-body ones. Since there are two different 
ways of factorizing, for consistency one must have 

Si(9i — 9 2 )Si + i(9i — 9 3 )Si(9 2 — 63) 
= S l+1 (9 2 - 93)3,(9, - 3 )5 l+1 (0! - 62). (3.2) 

The connection to the braid group is now obvious: 5(0) 
and 5(oo) obey the braid group relation l|2.1[) . In most 
known cases (and in the cases of interest here) 5(0) cx /, 
and S{-9)S{9) = I. Thus we have 

B = lim S(9) 

9^00 

B- 1 = lim S(-9) (3.3) 

8— »oo 

The matrix 5 = e lS ^e l9A Se lBB , where 8(6) is a func- 
tion of the rapidity 9, and A and B are diagonal 9- 
independent matrices. These factors arise in general to 
ensure that 5 has the correct properties under crossing 
symmetry and unitarity. Obviously, we need to remove 
the oscillating factors as 9 — > 00 to have a well defined 
limit. Both S(9) and the modified matrix S(9) satisfy 
the Yang-Baxter equation H3.2(l . 

The limit 9 — ► 00 in Eq. I|3.3fl also makes sense at an 
intuitive level. In order for the 5 matrix to be that of a 
loop gas, one should be in a limit where the mass m of 
the particles is small: otherwise, the loops would be high 
in energy and not dominate the partition function. When 



the particle mass is small, one can create particles of any 
rapidity 9i, and so the rapidity difference 9 = 9{ — 9 i+ i 
will typically be large. 

To conclude this section, we note that there are two im- 
portant additional steps to take in constructing a quan- 
tum loop gas having quasiparticles with non-Abelian 
statistics. The first is to find a Hamiltonian which has the 
ground-state wave function of Eq. . This can usually 
be done by a trick utilized by Rokhsar and Kivelsoni^. 
This trick is useful for any field theory with an explicit 
real action2244. For lattice models, there can be compli- 
cations, because one cannot always construct a Hamilto- 
nian which is ergodic in the Hilbert space. Nevertheless, 
in many cases of interest this procedure has been suc- 
cessful. 

The second additional step is to make sure that the 
excited states have braid relations which are those of the 
loops in the ground state. One way of doing this is to 
have a Hamiltonian so that the excited states are defects 
in the configuration space of loops. That is, a particle 
and an antiparticle over the ground state are connected 
by a strand. Thus when they are moved around each 
other, the non-locality due to the strand results in the 
braid relations described above. As is well known, this 
construction works in the Abelian case. 



IV. THE BRAID MATRICES 

In this section we give explicit expressions for the 
5 matrices and braid matrices associated with the 
Temperley-Lieb and BMW algebras described in section 
HH This will enable us in the next section to identify the 
two-dimensional classical field theories associated with 
these 1+1-dimensional quantum theories, so that we can 
construct quantum loop gases with the desired braiding. 

In the SU(2) case, the correspondence given in Eq. 
(|3.3|) means that we need to look for an 5 matrix which 
at infinite rapidities is of the form of Eq. (|2.3|) . Such an 
5 matrix has been known for quite some time^. It is 
straightforward to check that 

\e _ -\e 

Si{9) = 1 ~ g-V* - qe-» * ^ 

obeys the Yang-Baxter equation 1)3. 2f) for any value of the 
parameter A, as long as satisfies the Temperley-Lieb 
algebra, Ea. (|2.2|l . 

A number of related models have 5 matrices which can 
be written in the form of Eq. (|4.1() . The most famous is 
the sine-Gordon model. There are two different particles, 
the soliton (labeled +) and antisoliton (labeled — ), form- 
ing the spin-1/2 representation of the U q (sl2) symmetry 
of the model. The single-particle space of states Vi is two- 
dimensional, so that Ci (which acts on Vi®Vi + \) is a four- 
by-four matrix in this representation of the Temperley- 
Lieb algebra. Labeling the rows and columns in the order 



8 



++• H — i — I- , — gives 

/O 0\ 
s q 1 f . 

\0 0/ 

We have labeled this with a 6v because the the Boltz- 
mann weights of the six-vertex model can also be ex- 
pressed in terms of these e^. 

This particular representation e 6v of the Temperley- 
Lieb algebra, however, does not result in the braid ma- 
trices of the 2+1-dimensional theory, as it does not re- 
spect the Jones- Wenzl projector. Namely, consider k + 2 
strands in a row, i.e. the space V\ ® V% ® ■ ■ ■ ® Vk+2- 
Any k + 1 strands in a row must obey the Jones- Wenzl 
projector: we must restrict the space of states so that 
v ([k+i]/2) = v ([k+i]/2) = Q former acts non _trivially 

on the first k + 1 strands, the latter on the strands 
2 . . . k + 2). Now braid the k + 2th particle with the 
k+ 1st; the resulting configuration need no longer satisfy 
■p(l k + 1 }/ 2 ) _ q p or example, consider three strands in the 
k = 1 case, where configuration ( — | — ) in V\ V2 ® V3 
is part of the projected Hilbert space. If we braid the 
last two particles, the off-diagonal terms in e 6v result in 

a non-zero amplitude for the final state ( h). The 

latter state is projected out of the Hilbert space, since 
two S e = —1/2 states in a row are necessarily in a spin- 
1 representation. Thus the braiding does not commute 
with the projection: imposing the Jones- Wenzl projector 
violates unitarity. Obviously, we can't have this, so the 
only alternative is to conclude that ef" from Eq. (|4.2() 
cannot be used to build a braid matrix. 

Luckily, this issue is well-understood from a number of 
points of view. When k is an integer, there is another 
representation of the Temperley-Lieb algebra which pre- 
serves the projection. In the language of two-dimensional 
classical statistical mechanical lattice models, this repre- 
sentation is called the restricted solid-on-solid (RSOS) 
representation^^ 7 -. In the S matrix language, this rep- 
resentation describes the scattering of kinks in potential 
with k + 1 degenerate minima, as we discuss in Appendix 
lAl In the quantum-group picture, the RSOS represen- 
tation is obtained by "truncating" the states, so that all 
are in irreducible and indecomposable representations of 
U q {sl 2 ). 

The presence of the Jones- Wenzl projector means 
that we do not need to define the braid matrix on 
the full tensor product V® N , but only on the restric- 
tion/truncation/projection of the space V® N to the 

states obeying ■pj[ fc + 1 l/ 2 V® Ar = for all j. This re- 
stricted Hilbert space is our true space of states V(N). 
The states in V(N) are conveniently labeled in terms of a 
series of variables we call "dual spins" . (These variables 
are often called heights; we avoid this here to avoid confu- 
sion with the heights we discuss in the next section.) The 
dual spins take on integer values ranging from 1 . . . k+ 1, 
and live between the strands. Each strand is labeled by 



the two dual spins to the left and right of it, which must 
differ by ±1. The key effect of the restriction is the fact 
that dual spins only range from 1 to k + 1. This is a 
consequence of not allowing k + 1 consecutive strands to 
have spin (k + l)/2. 

A useful way of understanding the dual spins comes 
from treating each strand as being a spin-1/2 representa- 
tion of U q (sl2), and a dual spin r as being a spin-(r — 1)/2 
representation. Fix the first dual spin n to be have the 
value r% = 1, so it signifies the identity representation of 
U q {sl2). Our rules for dual spins mean that a strand next 
to 7"i separates this from a region of r-i =2; the region of 
dual spin 2 can be next to a region of dual spin rs = 1 
or r3 = 3, and so on. These are precisely the rules for 
taking tensor products in s^: we have ® 1/2 = 1/2, 
1/2 ® 1/2 = 0® 1, etc. In other words crossing a strand 
next to the dual spin r is like tensoring the spin-(r — 1)/2 
representation (the dual spin) with the spin-1/2 represen- 
tation (the strand). Thus we define the value (r— 1)/2 to 
be the total spin of the first i — 1 strands. Now imposing 
the Jones- Wenzl projector is easy. Forbidding the repre- 
sentation of spin (k+l)/2 is equivalent to forbidding the 
dual spin r = k + 2. Note that our earlier representation 
e 6v can be described in terms of dual spins as well: the 
+ particle increases the dual spin by 1 (moving left to 
right, say), while the — particle decreases it by 1. That 
the braiding from Ea. (|4.2|l can violate the Jones- Wenzl 
projection is obvious in the dual spin description: after 
braiding the value of the dual spin can reach k + 2 even 
if all the initial dual spins are below k + 2. 

To give the explicit RSOS representation of ej for k 
integer it is most convenient to label it by four dual spins 
r, s, t, u, each ranging between 1 and fc+ 1, with \r — s\ = 
\s — t\ = \t — u\ = \r — u\ = 1. The matrix elements of ti 
are then^S^ 




where [h] q = (q h — q~ h )/(q — q^ 1 )- The lines represent 
the strands; this picture represents what happens when 
the (sr) strand braids with the (ru) strand. After the 
intersection, the final state consists of the (si) and (tu) 
strands. Since the S matrix (and the corresponding braid 
matrix) are non-diagonal, the final state can be different, 
namely one can have r ^ t if s = u. A very important 
thing to note is that if s, r and u are between 1 and k+ 1, 
then the matrix elements for t = and t = k + 2 vanish 
because q k+2 = — 1. integer. Thus if an initial configu- 
ration satisfies the Jones- Wenzl projection, the final one 
does as well. 

Using the matrix of Ea.(|4.3p in Ea.l|4.ip gives the S 
matrix of an integrable 1+1-dimcnsional field theory; we 
identify this theory in the next section. Using the matrix 
of Ea. (|4.3|) in Ea. l|2.3|l (i.e. taking the 9 — > 00 limit of 
the S matrix) gives the braid matrix of particles associ- 
ated with SU(2)k Chern-Simons theory. This is also the 
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braiding of the quasiholes in the Read-Rezayi states in 
the fractional quantum Hall effect22iii4&. It is important 
to note that this representation, Ea. 14.3jl . only is useful 
for k an integer; otherwise the Jones- Wenzl projection 
cannot be satisfied (indeed, it does not exist). 

The space V® N for the (wrong) braid matrix Ij4.2|) has 
dimension 2 N . Since the actual Hilbert space V(N) is a 
subspace of V® N , its dimension must be smaller. Find- 
ing its size is a straightforward exercise done in many 
places; see e.g. Appendix^or Ref. Ell One finds for N 
large that the number of states grows as d . Thus the 
weighting per loop is indeed d, as the Temperley-Lieb 
algebra implies it should be. 

The results for 50(3) braiding and the BMW braiding 
can be derived from the Temperley-Lieb representation, 
Ea. (|4.3|) . There are two solutions to the Yang-Baxter 
equation whose 5 matrices turn into Ea. (|2.7|l in the in- 
finite rapidity limit 9 — > oo. In the classification of Ref. 
l49l one solution is associated with the fundamental rep- 
resentation of the U q {A 2 ) (A 2 is a twisted Kac-Moody 
algebra), while the other solution is associated with the 
spin-1 representation of U q (sl2). We will discuss in the 
next section how both solutions have been identified as 
5 matrices for two (related) field theories. The first so- 
lution is the most important for us. Written in terms of 
the generators X and E, it is 



Si 



q 2 e xe 



q -i e -\e 



qe 



xe 



q3 e \e 



-q-3 e -™ 3 



(4-4) 

Taking the 9 — > 00 limit yields the braiding matrix 
B so(3) in Eq .J27J. 

We can build a representation of the Ej and Xj from 
the ei of a Temperley-Lieb representation by using the 
relations of Eq. lj2.5jl . If the latter obeys the Jones- 
Wenzl projection, the resulting 50(3) BMW represen- 
tation does as well. In fact, we must use the represen- 
tation of Eq. H4.3fl . because it turns out that this yields 
the only unitary 5 matrix of the form of Eq. (|4.4I) . The 
representation of the states in terms of dual spins there- 
fore applies to the 50(3) case as well. However, here 
the rules for adjacent dual spins are different, since the 
strand is in the spin-1 representation of U q (sl 2 ). For 
dual spins l...k — 1, the rules are the same famil- 
iar ones from ordinary sli spin-1 representations: e.g. 
(h - l)/2 ® 1 = (h - 3)/2 © (h - l)/2 © (h + l)/2 for 
3 < h < k — 1. For dual spins k and k + 1 we have 
respectively (k ~ l)/2 1 = (k - 3)/2 © (fc - l)/2 and 
k/2 (g) 1 = (fc — 2)/2; in the latter, the representation of 
spin k/2 does not appear in the right-hand-side. Note 
that the states split into two subsectors, since even dual 
spins must always be adjacent to even dual spins, and 
odd are adjacent to odd. 

If one uses the rules given in Appendix 1X1 to count the 
number of states for N spin-1 strands, one finds it grows 
as (d 2 — 1) N at large N. Thus we can indeed interpret 
d 2 — 1 as the weight of an isolated loop, as the 50(3) 
BMW algebra implies. 



V. THE 1+1-DIMENSIONAL FIELD THEORIES 

In order to build our quantum loop gas, we need one 
more ingredient. This is to identify the underlying two- 
dimensional classical field theory, so that the wave func- 
tion of the 2+1-dimcnsional theory is given by Ea. lj3.ljl . 
We argued in Section ITTT1 that the 1+1-dimensional ver- 
sion of this underlying theory will have an 5 matrix 
whose 9 — > 00 limit gives the braid matrix. In this sec- 
tion we identify the field theories whose 5 matrices are 
those in the last section. 

To construct the quantum loop gas directly from the 
field theory, one needs to know the action of the two- 
dimensional classical field theory. As we will discuss in 
more detail below, for most of the theories of interest, 
the explicit action is fairly difficult to deal with. How- 
ever, as discussed in Appendix^ there are nice Landau- 
Ginzburg descriptions. Thus one can define a quantum 
Hamiltonian in this language, using the procedure dis- 
cussed in Refs. 44 and 129. 



A. The SU(2) case 

We s tart with the SU(2) case, arriving at the results 
of Refs. from a slightly different point of view. 

The 5 matrix, Ea. 14.lf) . with a in the representation of 
Eq. (|4.3|) , describes a field theory which can be defined in 
several different ways, which we describe here. 

One definition is as the continuum limit of an RSOS 
lattice models. The degrees of freedom of an RSOS 
lattice model are on the sites of a square lattice. The 
variables are called "heights" , and are integers ranging 
from 1 . . . k + 2. Heights on nearest-neighbor sites must 
differ by ±1. The Boltzmann weights for this model are 
those of regime III in Ref. ^(J This phase is ordered^. 
Each ordered state has only two heights present: one 
sublattice has all heights h while the other has all h + 1 . 
The excitations are the k different kinds of domain walls 
between the k + 1 different ordered states; each wall can 
be labeled by the two heights h,h±l it separates. 

It is quite simple to see qualitatively how the 5 ma- 
trix applies to this RSOS height model. The dual 
spins in the representation jUJl are identified with the 
ground states of the height model (there are k + 1 of 
each, with the rule that adjacent ones must differ by ±1). 
The strands are identified with the excitations of the lat- 
tice model, the domain walls. In the absence of defects 
the domain walls form non-intersecting loops, just like 
the particle worldlines we have described in detail above. 
The domain walls are indeed the objects whose scatter- 
ing is described by the 5 matrix. In the 1+1-dimcnsional 
picture, the excitations can be thought of as kinks, as dis- 
cussed in Appendix lAl 

Of course there are multiple lattice models with the 
same continuum 5 matrix. The RSOS height model has 
the advantage that it is integrable, and that the connec- 
tion of the 5 matrix to the lattice variables is quite intu- 
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itive. However, there is another lattice model in almost 
the same universality class. We say "almost" because 
some modifications are required if the space is a torus. 
This caveat does not affect the S matrix, and in any case 
we will not worry about the torus. The model with the 
same S matrix is called the lattice 0(n) loop model, be- 
cause at n integer it is 0(n) invariant. However, we are 
interested in the case n = d = q + q^ 1 = 2 cos(7r / (k + 2)) , 
so that |n| < 2. This model can be defined for all n as 
a gas of self- and mutually-avoiding loops on the honey- 
comb lattice with a weight n per loop, in addition to a 
weight per length of loop. By writing the S matrix for 
the 1+1-dimensional version of the 0(n) model in terms 
of generators obeying algebraic relations, one in fact can 
make at least formal sense of it for all values of n < 2, not 
just k integer—. These algebraic relations are equivalent 
to the Temperley-Lieb algebra™. The loops are inter- 
preted heuristically as the world lines of the particles. In 
these works, no explicit representation of the Temperley- 
Lieb algebra is necessary, but the Jones- Wenzl projection 
is required to obtain the correct answer for physical quan- 
tities on the cylinder—. Thus when k is an integer one 
can use the RSOS representation in the 0(n) model as 
well, although the physical interpretation of this in the 
lattice model is not very clear. 

(As a side remark, we note that the earlier representa- 
tion of Eq. H4.2f> does have a nice heuristic interpretation 
in the context of the 0(n) lattice model. One can formu- 
late this model as a model of oriented loops, where clock- 
wise loops get a weight q and counterclockwise loops get 
a weight q^ 1 . Despite the complex Boltzmann weights, 
the partition function remains real after summing over all 
orientations. The formulation in terms of oriented loops 
is useful because this can be mapped onto a model with 
local interactions, the six- vertex model with staggered 
Boltzmann weights 35 . The projection mentioned is nec- 
essary to get the correct weighting for loops which wrap 
around the cylinder. In this formulation, the + and — 
particles mentioned at the beginning of Section IIVI then 
correspond to the two orientations of the loop.) 

The description in terms of the OM lattice model 
is precisely that found in Refs. IM7U3CI The 2+1- 
dimensional lattice model discussed there has a ground- 
state wave function of the form of Ea. (|3.1|) . where the 
action S is precisely that of the O(n) loop model with 
n = d = q + q^ 1 . It was convincingly argued that this 
model indeed has fractional statistics, with braid matrix 
given by Ea. i|2.3f) . Thus by our S matrix line of argu- 
ment, we have arrived at the same conclusion. This is 
therefore strong evidence in favor of our conjecture in 
Section ITTT1 that the braid matrix of the 2+1-dimensional 
theory is related to the S matrix of the corresponding 
1+1-dimensional theory. 

The theories with the SU{2) RSOS S matrix can be 
formulated directly in the continuum, without need for 
the lattice descriptions given above. For general fc, how- 
ever, there is no simple field-theory action for these theo- 
ries, although a heuristic but useful Landau- Ginzburg de- 



scription is given in appendix^] They can also be defined 
in terms of constrained fermion models 5 ^ which realize 
the Goddard-Kent-Olive current algebra construction- 
It is difficult to obtain much information from this formu- 
lation, however. For our purposes, it is most convenient 
to define the field theories of interest as perturbations 
of a conformal field theory. One can define and indeed 
solve conformal field theories without a Lagrangian: the 
Hamiltonian and states are defined in terms of represen- 
tations of the Virasoro algebra. A massive field theory 
is defined by perturbing the conformal field theory by a 
relevant operator. As shown in Ref. l56l the S matrix of 
Eq. (|4.1() with a given by Eq. I|4.3[) is that of a perturba- 
tion of the conformal minimal model with central charge 



The desired S matrix describes the perturbation of the 
conformal field theory with p = k + 2 by its least relevant 
primary field (known usually as 'I'l.a), which has scaling 
dimension 2(p— l)/(p+ 1)). 

Before moving on to the SO(3) case, we wish to note 
another complication in the above picture. The first is 
that the S matrix of Ea. l|4.1[) applies to the 0(n) model 
in its dilute phase, where the energy per unit length is 
larger than the entropy, so that the loops cover a small 
part of the lattice. In order to get a purely topological 
field theory, the weight per unit length of loop must be 
1, so that no length scale is set for the loops— Such 
an 0(n) model for n < 2 is in its dense phase, where 
entropy wins and the loops cover a set of measure 1 of 
the lattice. However, the braid matrix is not related to 
the S matrix in the dense phase, but rather that of the 
dilute phase. (The S matrix in the dense phase has been 
studied, but due to the non-unitarity of the model, un- 
derstanding it precisely is a complicated and somewhat 
gruesome story.) The dense and dilute phases are not 
dual to each other; the former has algebraically decaying 
correlators, while the latter's decay exponentially. The 
same statements can be made in the context of the height 
models describing the perturbed minimal models. 

The way of understanding this complication is to re- 
member that in the dilute phase the arguments of Sec- 
tion IIIII suggest that S matrix really is describing the 
scattering of the excitations themselves, i.e. what hap- 
pens when two world lines braid. The braiding we are 
interested in is of the bare loops, not the renormalized 
excitations, and this is given by the S matrix in the di- 
lute phase. The S matrix in the dense phase is describing 
the excitations over the sea of dense loops, which is im- 
portant in the 1+1-dimensional case, but not of interest 
for the 2+1-dimensional braiding. The lesson is that the 
braiding should indeed be interpreted as that in the di- 
lute phase, even though the topological point is in the 
dense phase where loops proliferate. 
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B. The SO (3) theory 

Several field theories with the S matrix of Ea. (|4.4() 
were identified and discussed by Smirno\i 5 i The case 
we will focus on here corresponds to a perturbation of 
a minimal conformal field theories with central charge of 
Eq. (|5.1|) . However, for a given k, both the minimal model 
and the perturbation are different from the SU(2) case. 
This time, we have p = k + 1, and the perturbation is by 
the $21 operator. 

It is convenient to use the much-better-known interpre- 
tation of this field theory as the continuum description 
of the Q-state Potts model where Q is given byi&iS 

Q = d 2 = ( ? + <T 1 ) 2 = 4co8 a (r^Y 

The Potts model can defined for all Q in terms of its high- 
temperature expansion, where Q becomes a parameter. 
This definition does not have local Boltzmann weights 
for arbitrary Q, but for our special values with k integer, 
there is a lattice model with the same high-temperature 
expansion. This is found by using the original Temperley- 
Lieb result of writing the Potts transfer matrix in terms 
of generators obeying the algebra of Ea. l|2.2|l . and then 
using the RSOS representation, Eq.lj4.3p. of these gen- 
erators. At a particular coupling where the weights are 
isotropic, the lattice models are identical to the RSOS 
lattice models^ 6 ^ at their critical point. Thus the Potts 
critical point is also described by the conformal field the- 
ory with p = k + 1. Off the critical point, the S ma- 
trix of the Potts model with k integer is indeed of the 
form of Eq. H4.4|) 42 - 5, - 59 - 6 ° However, when Q is an integer 
(k = 2,4, oo), this 5* matrix is diagonal, so the braiding 
which follows from it is Abelian. For non-Abelian statis- 
tics, we need to use the Potts model for Q not an integer. 
We discuss this quantum loop gas in detail in Sect ion IvH 

As opposed to the 0(n) model for n ^ 1, the Potts 
model has a duality relating high to low temperature. 
On the lattice, this is a generalization of the Kramers- 
Wannier duality of the Ising modelSi* 6 ^. In the conformal 
field theory picture, there is a Z2 symmetry relating the 
perturbing operator $2,1 to — $2.1- All the operators 
appearing in the operator product expansion of $2,1 with 
itself are irrelevant, so perturbing by $2,1 and —$2,1 must 
be equivalent. The two signs of perturbation correspond 
to the low- and high-temperature phases, with the critical 
point being the self-dual point. 

This duality is a crucial ingredient in interpreting 
states in the 2+1-dimensional quantum model. We have 
stressed above how excitations with non-Abelian statis- 
tics can arise in quantum loop gases, where the ground 
state is a liquid state, i.e. a superposition of many states 
which does not break any symmetries. The key to under- 
standing how to do this here is to view what we referred 
to as the high- "temperature" phase in the classical statis- 
tical mechanical picture as a quantum disordered ground 
state which we picture to be a superposition of the exci- 
tations of the dual ordered phase, i.e. the excitations of 



the classical low- "temperature" phase. Recall that in the 
2+1-dimcnsional picture, the weights measure the am- 
plitude of a particular configuration in a wave function. 
This terminology contains somewhat of an abuse of lan- 
guage: the quantum system is not at high temperature, 
but rather the weights of the ground state corresponds 
to high temperature in the classical model. We are dis- 
cussing the properties of the quantum system only at zero 
(physical!) temperature. 

In the Abelian case, this can be seen quite clearly in 
Kitaev's model 5 . Here the underlying classical lattice 
model is the Ising model. This is therefore equivalent 
to both the SU{2)\ model (based on the O(l) loop gas) 
and our £0(3)2 model (the Q = 2-state Potts model). 
The loops are simply the domain walls between the Ising 
spins, which get a weig ht 1 = (V2) 2 - 1. The corre- 
sponding 2+1-dimensional model is topological when the 
action S in the wave function of Eq. I|3.1f) is of the Ising 
model at infinite temperature, where the Ising domain 
walls have zero energy per unit length and have prolif- 
erated. In the ordered phase, the order operators have 
expectation values, and the excitations are created by the 
disorder operators 2 *!. In the disordered phase, the disor- 
der operators get expectation values, and the excitations 
are created by the order operators. The lesson is that 
when there is a duality, the operator which creates exci- 
tations in one phase is the one which gets the expectation 
value in the dual phas o 63 i 64 . 

To conclude this section, we recall that there is another 
model which has an S matrix of the form of Ea. Q4.4[l . 
This is the tricritical Potts model, which in conformal 
field theory language corresponds to the minimal model 
with p = k + 2 perturbed by the $1,2 operator. One 
could presumably build quantum loop gases based on the 
tricritical models as well. Since in two dimensions, the 
tricritical point is unstable to perturbations toward the 
ordinary critical point, this would presumably hold as 
well in the 2+1-dimensional version. Thus such a quan- 
tum loop gas would be near a multicritical point as well. 

We also noted above that there is a second S matrix 
which reduces to _B S °( 3 ) in the 9 — > 00 limit. This S 
matrix is associated with a certain perturbation of the 
SU(2)k/U(l) "parafermion" conformal field theories 65 . 
(The perturbation is the Wi^i operator, where is the 
parafermion operators.) The physics is different for the 
two signs of this perturbation. For one sign, one obtains 
a massive phase, with this S matrix. For the other sign, 
one flows to the minimal model with the central charge 
given in Eq . (|5. 1|) with p = k + 1. This is precisely the 
critical point of the Potts models! Moreover, both critical 
points appear in the same RSOS lattice model (in the 
nomenclature of Ref. l46l the parafermion critical point 
separates regimes I and II, the minimal model separates 
regimes III and IV). Thus our interpretation is that this 
second S matrix is describing the same quantum loop gas 
in a region near another multicritical point. 
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VI. LATTICE MODELS 

We can combine all these ingredients to build quantum 
loop gases on suitable lattices whose excitations should 
have non-Abelian statistics. As discussed in Section ITTT1 
the strategy is to build a model whose ground state is 
given by a loop gas where the loops have the correct 
properties (e.g. a weight of d per loop in the SU(2) case, 
and a weight of d 2 — 1 = Q — 1 per isolated loop in the 
50(3) case). Such a lattice model for the SU(2) case 
was introduced by Freedman, Nayak and ShtengelStt. We 
repeat some of these arguments here, and then use the 
S matrix picture to define an analogous model for the 
SO(3) case. 

A. Criteria for the lattice models 

In all of the lattice models discussed here, the Hamil- 
tonian is of the Rokhsar-Kivelson form, meaning it can 
be written in terms of a sum of projection operators: 

H = ^\ i H i . (6.1) 

i 

The projection operators Hi = Hf are local but not nec- 
essarily commuting. The off-diagonal terms in the Hi 
must be ergodic, in the sense that any configuration can 
be mapped to any other (with the same values of any 
globally conserved charges) by repeated applications of 
the Hi. To obtain a desired ground state |\E'), one must 
find a set of operators Hi so that 

ffil*) = 

for all i. This means that the state |\&) is an eigenstate 
of H with energy 0. As long as all the coupling constants 
Xi are strictly positive, Aj > 0, this state |\&) is a ground 
state. We study models where the solution of this equa- 
tion can be written in the form of Eq. l|3.1|) : the basis 
elements of the Hilbert space can be thought of as a con- 
figuration in a classical two-dimensional lattice model, 
and the weight of this configuration can be expressed in 
terms of a local action. A key requirement that we will 
impose is that of locality i.e. that all the operators Hi 
act on a finite set of contiguous degrees of freedom. 

The degrees of freedom of the models in this section 
consist of a quantum two-state variable on each link of 
some two-dimensional lattice. We call these two states 
occupied and unoccupied. An occupied link corresponds 
to the presence of the strand, which one can think of 
as being in the spin 1/2 representation of U q (sl2) in the 
SU{2) case, and spin 1 in the SO (3) case. An empty 
site corresponds to the identity representation. What 
this means is that at each vertex, configurations appear- 
ing in the ground state must obey the corresponding 
fusion rules of £/ q (sZ 2 ) ■ 34 i 41 i 66 ' 67 i 68 ■ For example, three 
links, in states corresponding to representations r, s and 
t of U q (sl<i), touch each vertex of the honeycomb lattice. 



Configurations in the ground state must have the identity 
representation in the tensor product r ® s ®t. Thus in 
the SU(2) case, at each vertex must be touched by zero 
or two occupied links. In the 5*0(3) case, each vertex 
must be touched by zero, two, or three occupied links. 



1. The SU{2) Lattice Loop Models 

For the SU (2) case, one needs a set of Hi which annihi- 
lates states with the weighting rules of the 0(n) model. 
In other words, the ground state must consist of a su- 
perposition of configurations where the strands form self- 
arid mutually-avoiding loops. Moreover, each loop should 
have a weight d, and to be a purely topological ground 
state, there should be no weight per unit length. Pre- 
cisely, the criteria imposed on the configurations in the 
ground state at the purely topological SU(2) point are?2S 

1. The strands form closed non- intersecting loops: i.e. 
each vertex has either or 2 links with occupied 
links touching it. 

2. If two configurations are related by moving strands 
around, without cutting the strands or crossing any 
other strands, then these two configurations must 
have the same weight. In other words, two topologi- 
cally identical configurations have the same weight. 

3. If two configurations are identical except for one 
having a closed loop around a single plaquette {e.g. 
a loop of length 6 on the honeycomb lattice, length 
4 on the square), then the weight of the configu- 
ration without the single-plaquette loop is d times 
that of the one with it. 

The latter two properties are known as <i-isotopji Note 
that arbitrarily-sized loops are not directly required to 
have weight d. Rather, this property follows indirectly 
by combining the two latter properties: one can use cri- 
terion 2 to shrink a loop to its minimal size, and then use 
criterion 3 to remove it altogether while giving a relative 
weight d to the ground-state wave function. 

It is now straightforward to find the Hi annihilating a 
state with these properties by using locally-defined pro- 
jection operators. An explicit expression for the Hi in 
the SU{2) case on the honeycomb lattice can be found in 
Ref. l30i Since we will not need the explicit Hamiltonian, 
we will not give it here - it's rather ugly, but it does the 
job. In these models so far, d is a parameter that can 
take on any value, as the Jones- Wenzl projector is not 
imposed. To impose this projector on the ground state, 
one can add an energy penalty for configurations which 
violate the projection. This requires a fine-tuned inter- 
action involving a number of terms involving k + 1 spins 
(or strands) for level k. 
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2. The Lattice SO (3) Models 

To study the 5*0(3) theories for arbitrary k, we need to 
work harder. The appropriate lattice models are found 
by imposing criteria analogous to those of the SU (2) case, 
but adapted to the spin-1 loops. 

A typical loop configuration in the 5*0(3) model should 
look like that in Fig. |H1 The lines in this figure represent 



FIG. 8: A typical configuration in the spin-1 loop model 

"spin-1" particles, so that they correspond to the pro- 
jected double lines in the earlier Figs. Eland[7| We thus 
dub this the spin-1 loop model. We still require that the 
strands form closed loops. However, as opposed to the 
SU(2) case, we must now allow for trivalent vertices, i.e. 
the loops are now allowed to branch and merge. Thus the 
spin-1 loop model has branching loops, i.e. the configura- 
tions are nets^i. In the language of the quantum-group 
algebra U q (sl2), the reason for the trivalent vertices is 
that spin 1 appears in the tensor product of two spin-1 
representations. Equivalently, one can form an invariant 
from three spin-1 representations. Pictorially, this fol- 
lows from the presence of the 50(3) BMW generator X 
in Fig. This generator does not occur in the SU(2) 
case. In the quantum-group language, there are three 
generators here because the three representations of spin 
0, 1 and 2 appear in the tensor product of two spin-1 
representations. Precisely, the projection of two spin-1 
strands onto a spin-zero strand is E/(Q — 1), and onto a 
spin-1 strand is (X - E)/(Q - 2). 

In the spin-1/2 loop model for the lattice 0(n) 
modeffi, the Temperley-Lieb relation e? = dei implies 
that each loop receives a weight d. We can understand 
the somewhat more intricate analogous properties of the 
spin-1 loop model by using the 50(3) BMW algebra. 
The relation E? = (Q — l)Ej implies that isolated loops 
in the spin-1 model receive a weight of Q — 1 = d 2 — 1. 
Because trivalent vertices occur here, however, all loops 
need not be isolated. The projector onto spin-1 is pro- 
portional to X — E, so we associate this with two neigh- 
boring trivalent vertices, as indicated in Fig. Several 
properties of the loop gas follow from this. The relation 
(Xj — Ej) 2 = (Q — 2)(Xj — Ej) means that a configu- 
ration with a loop with just two lines emanating from 
it has a weight Q — 2 times the configuration with the 
loop removed. This is illustrated in Fig. ^| Moreover, 
because (Xj — Ej)Ej — 0, no graph can contain any loop 




= X- E 



FIG. 9: Two trivalent vertices in the SO (3) BMW algebra 




(Q-2) 



FIG. 10: Removing one loop in the spin-1 loop model 



with just one external line attached to it. We call such a 
forbidden loop a "tadpole" . 

We must work harder to find the weight of more com- 
plicated configurations in the loop model. To make the 
answer precise, we use a two-dimensional classical lattice 
model which has a loop expansion with the desired prop- 
erties. As we discussed in Section^ for the spin-1 loop 
model this should be the Q-state Potts model, since its S 
matrix gives the desired braid matrix. The desired loop 
expansion is the Zow-temperature expansion of the Potts 
model. 

Let us first describe the low-temperature loop expan- 
sion for Q integer, where the Potts models are defined by 
placing a "spin" cti taking values 1 ... Q at the sites i of 
a lattice. As it is well known, the interaction for a Potts 
model depends only on whether nearest-neighbor spins 
are the same or different, so that the Boltzmann weight 
for a link with spins <Xj and (jj at its ends is 

The low-temperature expansion is given by first express- 
ing each configuration of Potts spins in terms of its do- 
main walls. The domain walls reside on the links of the 
dual lattice, and separate regions on the direct lattice 
of spins of different values. Each link crossing a domain 
wall has weight e~ K , while each link without a domain 
wall has weight 1. The Boltzmann weight of a given spin 
configuration depends only on the length L of its domain 
walls: L is the number of links on the dual lattice with 
walls on them. The Boltzmann weight of a configuration 
is then e~ KL . A weight of 1 per unit length corresponds 
to infinite temperature in this classical lattice model. 

By definition, the domain walls must form loops sur- 
rounding groups of like Potts spins. These loops can 
intersect, but no tadpoles can occur. For example, the 
trivalent vertex given in Fig. ^2 occurs for Q > 2. 

Different configurations of spins can have the same 
domain-wall configuration: e.g. there are Q configura- 
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cates that the corresponding two loops are adjacent, and 
a dashed line indicates that two formerly independent 
loops are now merged (i.e. the occupied links separating 
them are removed). This is fairly obvious in the coloring 



FIG. 11: A trivalent vertex in the 3-state Potts model 



FIG. 12: The recursion relation for the chromatic polynomial 



tions with no domain walls, and Q(Q — 1) configurations 
with a loop of length 4 surrounding a given site. In gen- 
eral, the number of spin configurations which have the 
same loop configuration C is called the number of Q- 
colorings Xq(£) ( see e -9- Ref . 16^) . Imagine each region 
of like spins to be shaded some color. The number \q 
is then the number of ways this shading can be done 
with Q colors so that no two adjacent regions have the 
same color (regions which meet only at a point are not 
considered to be adjacent). 

The partition function of the Potts model can therefore 
be written as 

c 

where the sum is over all distinct loop configurations: the 
multiple spin configurations with the same loop config- 
uration is accounted for by the factor Xq(£)- A typical 
loop configuration C looks like that in fig. [5] Because 
Xq{£) vanishes for any configuration with a tadpole, or 
a strands with dangling end, we need not include such 
configurations in the sum. This expansion is a useful de- 
scription of the ferromagnetic (K > 0) Potts model at 
low temperature. It is important to note that this is not 
the only loop expansion of the Potts model: another ex- 
pansion is in terms of the (self- and mutually-avoiding) 
loops surrounding the clusters in the high-temperature 
expansioni^iiiiiLZil. 

The low-temperature expansion of the partition func- 
tion of the Potts model, Eq.JfOJ), applies to any Q when 
Xq{£) is the chromatic polynomial of the graph dual to 
£iS The graph dual to C is defined with a node cor- 
responding to each loop, and a line between two nodes 
when the corresponding loops share a boundary. In terms 
of the Potts spins, each node in this graph corresponds to 
a region of like spins, and a line between two nodes means 
that corresponding regions are adjacent. The chromatic 
polynomial reduces to the number of colorings of the 
graph when Q is an integer, but can be defined for all 
Q by a recursion relation. Consider two nodes connected 
by a line I (i.e. two loops sharing a boundary in the orig- 
inal picture). Then define T>iC to be the graph with the 
line deleted, and CiC to be the graph with the two nodes 
joined into one. Then we have 

X Q (C)=x Q (V l C)~x Q (C l C). (6.3) 

We represent this pictorially in fig. 1121 where a node rep- 
resents each loop, a solid line between two nodes indi- 



description: Xq (Pity includes all graphs in Xq(£), but 
also has graphs where the two nodes connected by line I 
have the same color. There are Xq(C;£) of the latter so 
we need to subtract these off to get the recursion relation. 
For any loop configuration £, one can apply l|b.3fl repeat- 
edly until one reaches graphs with all isolated nodes. A 
graph with Af isolated nodes has Xq = ■ We will give 
explicit examples of how this works in Subsection IVIBI 
The criteria for the Potts loop model to describe the 
ground state of the £0(3) loop gas on the lattice are 
therefore 

1'. The strands form closed loops, but now we allow 
trivalent vertices. 

2'. If two configurations are related by moving strands 
around without cutting the strands or crossing any 
other strands, then these two configurations must 
have the same weight. In other words, two topologi- 
cally identical configurations have the same weight. 

3'. Each loop configuration C receives a weight Xq{£)- 
For example, if two configurations are identical ex- 
cept for one having a closed loop around a single 
plaquette (a loop of length 6 on the honeycomb 
lattice, length 4 on the square), then the weight of 
the configuration without the single-plaquette loop 
is Q — 1 times that of the one with it. 

Criterion 2' is the same as criterion 2 in the SU (2) case; 
this is the requirement of topological invariance. Crite- 
rion 1' is the generalization of criterion 1, allowing for 
trivalent vertices in the SO(3) case. Criteria 3' is the ap- 
propriate generalization of criterion 3. However, imple- 
menting this using a local Hamiltonian requires a little 
work, which we will now describe. 

B. A Hamiltonian yielding SO(3) statistics 

In the previous Subsection we set out the criteria which 
the ground-state wave function for the SU(2) and SO(3) 
models must obey. Here we describe a Hamiltonian of the 
form of Ea. (|6.1(l for the SO(3) case, which has a ground 
state with the weights of the (low-temperature) Potts 
loop gas. This is tantamount to finding a set of projection 
operators Hi which annihilate the desired ground state, 
and which result in an ergodic Hamiltonian. 

There are two types of Hi operators in our Hamilto- 
nian. The simplest type have purely potential terms, 
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diagonal terms where Hi = 1 on some basis elements of 
the Hilbert space, and zero on the remaining elements. 
Such terms thus allow us to satisfy criterion 1': we give a 
positive potential to any vertex which has only one occu- 
pied link touching it. Since the ground state has energy 
zero by construction, any state on which these Hi are 
non-zero cannot be part of the ground state, as long as 
there are no off-diagonal terms which mix this state with 
an allowed one. 

The second type of term contains off-diagonal ele- 
ments, and are needed to ensure that different basis ele- 
ments have the desired relative weighting in the ground 
state. A state of the form |£) + a\C) is annihilated by 

H = -ix) (6.4) 

so that e.g. Hi\C) = a\C) — \C). Since we know each 
state £ in the ground state receives a weight Xq(£), we 

need a = Xq{£) / Xq{£) ■ 

However, if we include an element Hi in the Hamil- 
tonian for any pair of states £, £', the Hamiltonian will 
be non-local for several reasons. An obvious one is that 
if we include an Hi for any pair of configurations, the 
off-diagonal terms are clearly non-local. We can easily 
solve this problem by setting Aj = for any Hi involving 
an £ and a £' whose differences arc non-local. In other 
words, we only allow off-diagonal terms in H which map 
a given £ to an £' which differs from £ only in some 
small neighborhood (say the links on a given plaquette). 
While this is a necessary condition for a local Hamilto- 
nian, it is not sufficient for this model. The reason is that 
evaluating Xq(£) for a given £ is a non-local operation: 
it requires knowing the entire cluster of occupied links. 
However, even though the overall Xq(£) needs to be de- 
termined globally, the ratio Xq(£)/xq(£') m some cases 
depends only on local difference between £ and £', not 
their global form. Since the Hamiltonian only depends 
on this ratio, we can find a local Hamiltonian if we can 
find such pairs £ and £'. 

Let us first describe how to implement criterion 2', 
which says that if two configurations £ and £' are topo- 
logically identical, then they have the same weight in the 
ground state. By definition, if £ is topologically identical 
to £', we have Xq(£) = Xq(£'). Thus including Hi of 
the form of Eci. (|6.4[l . with a — 1 will insure the proper 
weighting. These terms will be local if we require £ and 
£' to not only be topologically identical, but completely 
identical except on the links around one plaquette. 

To make these terms in the Hamiltonian more specific, 
let us work henceforth on the honeycomb lattice, so that 
we do not have to worry about loops which touch at only 
a point. Consider a single plaquette, where some but 
not all of its six links are occupied. The simplest pos- 
sibility allowed by criterion 1' is then for two of the six 
outside links touching the plaquette to be occupied, as 
in all the configurations in Fig. El For each configura- 
tion of the two outside occupied links (15 possibilities in 




FIG. 13: Plaquettes with two external occupied links 



all), there are two topologically-identical configurations 
on each plaquette. We thus include in H the a = 1 pro- 
jectors which include flips between the two topologically- 
identical configurations; two of the flips are illustrated in 
Fig. El These Hi are local, involving only states on 12 
links: the six on the plaquette and the six touching it. 

This idea can readily be generalized to plaquettes with 
more of the outside links occupied. If there are three 
outside links occupied, then there are three topologically- 
identical configurations on each plaquette, as illustrated 
in one case in Fig. El We thus include a = 1 projec- 




FIG. 14: Plaquettes with three external occupied links 



tors which flip between any pair of topologically-identical 
configurations. For four outside lines, there are two pos- 
sibilities. The first type of configuration shown in Fig. 1151 
has no topologically-identical partner, while the second 
type has one. We thus include a — 1 projectors for all 




FIG. 15: Plaquettes with four external occupied links 



configurations of the latter type. For five or six outside 
lines, we include no projectors. By repeatedly applying 
the Hi described in Figs. El El and El we implement 
criterion 2'. 

As if this Hamiltonian weren't already complicated 
enough, we now need to implement criterion 3'. The 
Hi described in Figs. El & U map between topologically- 
identical configurations. To map between topologically- 
distinct configurations, we need still more Hi operators. 
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These are still of the form of Eg. l|6.4[l . but to ensure dif- 
ferent configurations have the correct relative weight, the 
a are not necessarily equal to 1. Again, we focus on a 
single plaquette, but here we consider only plaquettes 
with all six links occupied. The Hi depend of course on 
which of the outside links are occupied. The cases with 
0, 1, 2, and 3 occupied outside links and all internal links 
occupied are easy to implement. We have: 

0. When there are no outside links occupied, then a pla- 
quette with all six links occupied forms an isolated 
loop. If we remove this loop, the resulting configu- 
ration has weight Q — 1 relative to the configuration 
with the loop. This can be implemented with an 
Hi with a = Q — 1, where C is the configuration 
with the isolated loop and £ is the configuration 
without it. This is illustrated in fig. 1161 

FIG. 16: Removing or adding an isolated loop 



1. If there is just one occupied link, this is a tadpole, 

and is forbidden in the Potts loop expansion and 
hence the ground state. Thus we add a potential- 
only term (i.e. H = 1 on a tadpole). Tadpoles 
comprised of loops larger than a plaquette end up 
being forbidden by using the isotopy: applying the 
Hamiltonian enough will shrink a given loop to a 
single plaquette, and potential here will then ex- 
clude it from the ground state. 

2. If there are two occupied outside links connected to 

loop, the configuration with the loop removed has 
relative weight Q — 2, as illustrated earlier in Fig. 
ITU1 This is implemented by Hi with a = Q — 2, 
where C is the configuration with the loop, and 
C is a configuration with the same external lines 
but without the complete loop. As illustrated in 
Fig. 1131 for a given pair of occupied outside links, 
there are two allowed configurations on the plaque- 
tte with no loop. We can include a Hi(a = Q — 2) 
for either or both of these two allowed configura- 
tions. 

3. If there are three occupied outside links connected to 

loop, the configuration with the loop removed has 
relative weight Q — 3. We therefore use an Hi with 
a = Q — 3. Here £ can be any one of the three 
configurations illustrated in Fig. 1141 

A Hamiltonian comprised of the Hi we have con- 
structed so far has ground states with the correct relative 
weightings. However, it is not ergodic: there are multiple 
distinct ground states which are not related by any of the 



above off-diagonal terms. For example, a configuration 
with four adjacent plaquettes with all six links occupied, 
as illustrated in Fig. El below, is annihilated by all the 
Hi we have discussed so far. Thus we need more terms 
in H so that this state alone is not a ground state. Re- 
moving one of these loops (all of which have at least four 
occupied outside links) is not as simple as with three or 
less occupied links. We need to use the recursion relation 
Eq. (|6.3(l for the chromatic polynomial to find Hi which 
remove these loops. 

To find these terms, it is convenient to use the graph- 
ical representation of a loop configuration, as defined 
above. Knowing the graph of C is sufficient to find its 
chromatic polynomial Xq(£)- In terms of these graphs, 
the recursion relation Ea. (|6.3|) can be represented as in 
Fig. H21 above. We can use this relation to easily rederive 
the Hi acting on plaquettes with 0, 2 and 3 occupied out- 
side links. For example, the graph for an isolated node is 
precisely that on the left-hand-side of Fig. Applying 
Ea. (|6.3|) once, and then using the fact that an isolated 
node gives a factor Q to the chromatic polynomial, gives 
the desired relative weighting Q — 1. The equalities are 
meant as between the corresponding chromatic polyno- 
mials. 

Finding Hi which remove a loop with four external 
lines is trickier. One can apply the recursion relation, 
but one can get graphs which do not correspond to any 
configuration on the honeycomb lattice given by chang- 
ing links on the plaquette from occupied to unoccupied. 
Consider the first graphical representation of four adja- 
cent hexagons in Fig. 1171 In this graph, we allow the 




FIG. 17: Graphical representations of two configurations with 
five loops 



nodes B, C, D and E to be attached to other nodes, but 
node A touches only the four in the picture. We apply 
the recursion relation once to remove one of the lines at- 
tached to node A. This gives a perfectly valid relation 
for the corresponding chromatic polynomials, but there 
is no loop configuration corresponding to the graph with 



17 



one line removed. For example, say we remove the line 
from A to D: no loop configuration corresponding to this 
graph can be drawn on the honeycomb lattice. 

To define a Hamiltonian, we need a relation between 
valid loop configurations, not just different chromatic 
polynomials. We can relate the two loop configurations 
in Fig. El I n the first graph, we use the recursion rela- 
tion to remove the line from A to D and then the line 
from A to E. Now node A has lines only to nodes B 
and C, and corresponds to the situation of Fig. El We 
can now remove node A altogether, multiplying the re- 
sult (the square involving B,C,D and E) by Q — 2. This 
square defines a chromatic polynomial, but there is no 
corresponding loop configuration. We can however, re- 
late it to the second configuration in Fig. El the same 
square graph arises from using the recursion relation to 
remove the link from C to E. Combining the two, we 
obtain the relation in Fig. El with a = Q — 2 and b = — 1. 
Note that the first and last configurations are topologi- 




FIG. 18: Removing one loop of a four-hexagon cluster 



cally identical, so we can choose any value of a if we set 
b = Q — 3 — a. Thus a more symmetric Hamiltonian will 
result if we choose a = b = (Q — 3)/2. As a check on 
this relation, we can look at the special case where none 
of the nodes A, C, D, E are connected to any others (i.e. 
the four hexagons A, C, D, E are surrounded by region 
B). One can then easily verify that both sides are equal 
to (Q-l)(Q-2)(Q-3) 2 . 

We can then include an Hi with £ the sum of loop 
configurations on the right-hand-side of Fig. El One 
can define analogous relations to define new TLi which 
further reduce the number of loops on the right-hand- 
side of Fig. El Although we have not proven so, we 
believe that proceeding in this fashion one can define 
local Hi which result in an ergodic Hamiltonian. This 
should also be possible on other lattices, but presumably 
be even more complicated, since one must worry about 
loops which touch at only a point. Obviously, there is no 
conceivable way such finely-tuned Hamiltonians could be 
realized in nature. However, the fact that they are local 
makes it at least possible that there exists a more natural 
Hamiltonian in the same universality class. 



VII. THE PHASE TRANSITION 

Since we have identified the classical field theories 
which describe the ground state of the quantum loop 
gas, we can determine the phase structure of the latter. 
There is a major subtlety in doing so. So far, we have 
been discussing properties of the wave function tp. How- 
ever, in determining equal-time correlation functions of 
the quantum system, the functional integral is weighted 
by \ip\ 2 . In terms of the action of the corresponding clas- 
sical model, we have from Ea. (|3.1|) . 

|^( s )| 2 =e - s W- s *W 

for a configuration s in the classical model. Thus when 
we are computing for example correlators in the quan- 
tum model, we need to square the Boltzmann weights of 
the classical model. In a model where the loops them- 
selves are the degrees of freedom (such as those discussed 
in Section IVT1 the ones without the Jones- Wenzl projec- 
tor imposed), this means the weight per loop must be 
squared. In particular, this means loops in the SU (2) 
model get weight d 2 , while isolated loops in the 50(3) 
model have weight (d 2 — l) 2 . 

The phase structure depends on the weight per length 
of the loops. Let us first review what happens in a quan- 
tum eight-vertex model*** 2 ^. There the degrees of freedom 
in the ground state are arrows on the links of the square 
lattice obeying the eight-vertex condition (an even num- 
ber of arrows pointing in at each vertex) . The loops are 
given by following (say the up and right-pointing) arrows 
around the lattice; the eight-vertex condition ensures the 
loops are closed. One can rewrite these degrees of free- 
dom as Ising spins at the center of each plaquette, and 
the loops then form domain walls around the spins. The 
purely topological (Kitaev) point 5 , corresponding to our 
SU(2)x and 0(3)2 models, has equal amplitude for all 
configurations. This is indeed infinite temperature for 
the Ising spins. By including a weight per length of the 
loop, one can move away from the Kitaev point. As de- 
tailed in Ref . |2^, there is eventually a phase transition to 
an ordered phase. 

The interesting question to answer is if this is a critical 
point or a first-order phase transition. In the quantum 
eight-vertex model of Ref. |2^, the transition is second 
order. In the analogous transition in the SU(2) case&S, 
each loop will get a weight d 2 , so we can view this as 
a loop model with n c g = d 2 . The 0(n c ff) model has a 
critical point only when n e ff = 4cos 2 (7r/(fc + 2)) < 2, so 
this occurs only for k = 1 and k = 2. Only the latter is 
non-Abelian. 

The analogous result in the classical Potts model is 
that the phase transition at the self-dual point is second- 
order if Q < 4, first-order for Q > 4. This result, how- 
ever, cannot instantly be applied to the 2+1-dimcnsional 
case, because of the weighting by \ip\ . the phase struc- 
ture is that of the classical loop gas where each configu- 
ration is weighted by (xq(C)) 2 . This loop gas does not 
seem to have been studied before, so we do not know the 
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answer. However, we can make a simple conjecture. The 
configurations of this squared loop gas are of course the 
same as those of the Potts loop gas; only the weight per 
loop has been squared. One might therefore hope that 
the phase transition of the squared loop gas is in the same 
universality class as the Potts loop gas at some Q c ff ■ The 
simplest possibility is to assume that Q c f[ gives isolated 
loops in the squared loop gas the correct weight (Q — l) 2 . 
This amounts to a <3 c ff of 

Q cS -i = (Q-i) 2 = (d 2 -i) 2 . (7.1) 

The weight per unit length is also changed, but this only 
changes the location of the critical point, not its type. 
Thus we conjecture that there will be a critical point 
in the quantum loop gas if Q c s < 4, and a first-order 
transition otherwise. 

If this conjecture is true, critical points occur for 
k = 2,3 in the 50(3) model (as noted above, k = 1 is a 
trivial theory here). The k — 3 model has non-Abelian 
statistics. As noted at the end of section Hi CI this "Lee- 
Yang" model is the simplest model of non-Abelian statis- 
tics, since it has only one kind of strand. We have thus 
shown that in the quantum loop gas which realizes par- 
ticles with these statistics, there exists a critical point 
separating the topological phase from an ordered one. 

However, there's a catch. The Jones- Wenzl projector 
needs to be imposed separately to the loop models dis- 
cussed in Section lvTl at least when space is an annulus or 
a torus. This presumably amounts to a relevant operator 
at the critical poinliSi. Thus in the non-Abelian case even 
at k = 3, reaching the critical point from the topological 
phase requires tuning another parameter away. 

This is in harmony with some of our early observations. 
At the end of Section Til Bl we noted that for the 50(3) 
model with k — 3, one can use the Jones- Wenzl projector 
to remove X vertex, leaving only self-avoiding loops. The 
50(3)3 model with the projection is therefore equivalent 
(at least locally) to the 5/7(2)3 model. The latter does 
not have a critical point, even without the Jones- Wenzl 
projector, because k = 3 in SU (2) case corresponds to an 
0(n e s) loop model with n c g > 2. This is a strong indica- 
tion that imposing the Jones- Wenzl results in a relevant 
perturbation of the critical theory 
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APPENDIX A: THE LANDAU-GINZBURG 
DESCRIPTION OF THE 1+1-DIMENSIONAL 
THEORIES 

In this appendix we explain how the field theories dis- 
cussed in section [V] and their 5 matrices discussed in 
section llVl have a nice description in terms of a Landau- 
Ginzburg effective field theory. 

As shown in Ref. ITU a simple Landau-Ginzburg de- 
scription of the minimal model of conformal field theory 
with central charge of Ea. (|5.1() is in terms of a single 
scalar field 4> with potential ip 2 ^ p ^ 1 ' > . The critical point 
of the Ising model is the p = 3 case, and the tricritical 
point is the p = 4 case, and the corresponding <f) and (j) 6 
potentials have long been known. The critical point of 
the SU(2)k case (the continuum limit of the 0(n) or re- 
stricted height models) has p = k + 2. The critical point 
of the SO(3)k case (the continuum limit of the Potts or 
dilute Aj. + i models) has p = k + 1. 

The primary fields of the conformal field theory also 
can be described in terms of 4>, so the massive field the- 
ories of interest have a Landau-Ginzburg description as 
well. For the SU{2)k case, the effective description of 
the $1,3 operator is (j) 2 ( p ~ 2 \ For the Ising model, this is 
indeed the usual <fi 2 mass term. As is well known, this 
field theory has a single quasiparticle with 5 = — 1. In 
this case, must act on a one-dimensional space, and 
it is just a number; for p = 3, e$ = 1. This also follows 
from imposing the Jones- Wenzl projector: for k = 1, one 

imposes V^ 1 \d = 1) = I — = 0, which indeed gives 
ei = I. Thus the statistics for SU (2)i is AbelianSi 

For higher values of p, perturbing a 2 ( p_1 ) potential 
by cf) 2 ( p ~ 2 ' seems to induce a flow to the minimal model 
with p decreased by 1. This flow is indeed known to oc- 
cur for perturbations of one sign of $1,3, see Ref. |73 
In the 0{n) language, this perturbation corresponds to 
flowing into the dense phase. In the models of Ref. 
this is called regime IV. This is not what we want; this 
is a massless field theory with the problematic 5 ma- 
trix. Instead, we want to perturb by the same opera- 
tor $1.3 but with the opposite sign; the perturbations 
are not the same because there is no symmetry in the 
conformal field theory (except in the Ising case) which 
sends $1,3 to —$1,3. With this sign of the perturbation, 
the potential must renormalize to include extra terms. 
Since we know from the exact results^S^S that the lattice 
models have p — 1 ground states, we must fine-tune the 
potential to achieve this. This is very familiar from the 
tricritical Ising model p = 4. One sign of the perturba- 
tion moves the system along its first-order transition line, 
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while the other sign causes a flow to the ordinary Ising 
critical point. The Landau-Ginzburg potential along the 
first-order line indeed consists of a (j) 6 potential tuned to 
have three degenerate minima. 

Let us focus on p — 4 case in more detail. The Landau- 
Ginzburg potential for the tricritical Ising model along its 
first-order transition line along is <ft 6 + acjr + b<j) 2 , with 
b = —6a 2 so that there are three degenerate minima at 
= 0, ±\/2a- With such a potential, the low-energy con- 
figurations in the two-dimensional classical model consist 
of regions of these three vacua. The loops are domain 
walls between different vacua. When there are restric- 
tions on which vacua can be adjacent to each other, the 
allowed domain walls are restricted as well. What hap- 
pens here is that the vacuum +v / 2a is not allowed to 
be next to the vacuum — \[2a. Hence, there must be a 
region of vacuum in between. In the 1+1-dimensional 
description in terms of quasiparticles, the particles are 
kinks interpolating between adjacent vacua, and the do- 
main walls are their world lines. 

One can now count the "number" of particles. Say 
the left end of the system is in the vacuum 0. Then the 
space of states for one particle V(l) is two-dimensional; 
it consists of a kink going from to +V2a, and one go- 
ing from to — \/2a. However, the space V{2) is also 
two-dimensional: the reason is this restriction that the 
two vacua ±\/2a cannot be next to each other. Thus, 
if we have a kink going from to +^/2a, the next kink 
can only go back to again. The dimension of the al- 
lowed Hilbert space for N "particles", V(N) (with the 
boundary condition of vacuum at one end) , is therefore 
2[( Ar + 1 )/ 2 l ; where [x] is the integer part of x.This reduc- 
tion of the Hilbert space, with respect to the standard 
free particle Fock space, is the hallmark of non-Abelian 
statistics. Using the standard definition of the "number 
of particles" as ln(dim(V(A)))/7V, here we find \/2 This 
is indeed the correct value of d for k = 2, and gives a 
precise meaning to the statement that d, the weight per 



loop, can be interpreted as counting the "number of par- 
ticles going around a given loop" (albeit this "number" 
is y/2 !). 

Thus we see how these restricted kinks in 1+1 dimen- 
sions go hand-in-hand with non-Abelian statistics in 2+1 
dimensions. In the Moore- Read theory of the v = 5/2 
fractional quantum Hall effect, the quasiparticles have 
non-Abelian statistics, and the number of states obeys 
the same formula 2^ N / 2 \ where N here is the number of 
quasiholes ( "non-abelions" 

For general SU(2)k theories, we consider a 2 ( fc+1 ) po- 
tential with k + 1 minima. (The explicit potential can 
be written out in terms of Chebyshev polynomials if de- 
sired.) The kinks interpolate between adjacent vacua, 
so that their world lines are domain walls in the two- 
dimensional picture. To count the "number" of such kink 
configurations for large N, one applies the same proce- 
dure as above, and obtains 2cos N (ir/(k + 2)). In the 
Read-Rezayi generalizations of the Moore-Read theory 
of the fractional quantum Hall effect, the statistics and 
the number of states is the same^i. 

For the 0(3)k case, the perturbation of the conformal 
minimal model is different. The $21 operator is </> p_1 , but 
the resulting Landau-Ginzburg description is not as use- 
ful. However, in section ^ we showed how to describe 
representations of the 50(3) BMW algebra by fusing 
together representations of the Temperley-Lieb algebra, 
and the analog is possible here. We can describe the 
kinks in the Potts models as bound states of the kinks of 
the 0(n) model, although we believe that in the S matrix 
context this is a formal device without physical signifi- 
cance. Anyway, we consider the same potential with fc+1 
degenerate minima, but the 0(3)^ kinks are comprised 
of two of the SU(2)k kinks. More precisely, we consider 
all configurations made up of two SU(2)k kinks bound 
together, and then subtract off the identity. This proce- 
dure yields the correct fusion rules for spin-1 particles, 
and the correct "number" of kinks d 2 — 1. 
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